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We examine the validity of the harmonic approximation, where the radio-frequency ion trap is 
treated as a harmonic trap, in the problem regarding the controlled collision of a trapped atom and 
a single trapped ion. This is equivalent to studying the effect of the micromotion since this motion 
must be neglected for the trapped ion to be considered as a harmonic oscillator. By applying the 
transformation of Cook and Shankland we find that the micromotion can be represented by two 
periodically oscillating operators. In order to investigate the effect of the micromotion on the dy- 
namics of a trapped atom-ion system, we calculate (i) the coupling strengths of the micromotion 
operators by numerical integration and (ii) the quasienergies of the system by applying the Floquet 
formalism, a useful framework for studying periodic systems. It turns out that the micromotion is 
not negligible when the distance between the atom and the ion traps is shorter than a character- 
istic distance. Within this range the energy diagram of the system changes remarkably when the 
micromotion is taken into account, which leads to undesirable consequences for applications that 
are based on an adiabatic process of the trapped atom-ion system. We suggest a simple scheme 
for bypassing the micromotion effect in order to successfully implement a quantum controlled phase 
gate proposed previously and create an atom-ion macromolecule. The methods presented here are 
not restricted to trapped atom-ion systems and can be readily applied to studying the micromotion 
effect in any system involving a single trapped ion. 
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I. INTRODUCTION 

The dynamics of a trapped atom-ion system has re- 
cently been studied in great detail [1, i2j , and it shows an 
interesting application in realizing a two-qubit gate for 
quantum computing |3j. The system is composed of a 
single atom trapped in a harmonic trap interacting with 
an ion trapped in a radio-frequency (rf) trap. 

In these studies, the rapid motion of the ion on a short 
time-scale — the micromotion — is averaged out. This 
procedure is referred to in the literature as the harmonic 
approximation since it produces an effective harmonic 
motion of the ion, as if the ion were trapped in a time- 
independent harmonic trap. In fact, many of the propos- 
als for applications of this system are derived based on 
the harmonic approximation. However, there is a concern 
about the validity of this approximation since it is known 
that the kinetic energy of the micromotion is compara- 
ble to that of the ion's harmonic motion [i] . Motivated 
by such concerns, we present here two approaches, as 
mentioned in the Abstract, to studying the effect of the 
micromotion in trapped atom-ion systems. 

For this purpose we make use of the Floquet formal- 
ism [S]. The potential of the rf electric field used to 
trap ions depends periodically on time. When the po- 
tential is periodic, Floquet theory provides a powerful 
tool for treating the exact dynamics of the quantum sys- 
tem. It enables one to compute the so-called quasiener- 
gies and quasienergy states, also referred to as Floquet 
states, which are the analogs of the eigenenergies and 
eigenstates of a time-independent system. Hence, it of- 
fers a quantum-mechanical treatment of the micromotion 
problem in any system involving a single trapped ion. 
By studying the exact quasienergies and Floquet states 



we obtain valuable information about the role of the mi- 
cromotion in such systems. Its effect will be then de- 
duced by comparing the exact dynamics, as derived by 
the Floquet formalism, with the approximate one given 
by the harmonic approximation. In particular, this com- 
parison reflects on the validity of various proposals for 
applications using controlled collisions of trapped atoms 
and ions. Although we demonstrate the method for the 
system of interacting trapped atoms and ions, it can be 
readily applied to any system in which the trapped ion 
is coupled to an external time-independent subsystem. 

The article is organized as follows. In Sees. [II] and 
|III| we describe the ion's micromotion and the trapped 
atom-ion systems. Section IV is on the Floquet formal- 



ism and how we use it to study the effect of a period- 
ically oscillating potential on the energy structure of a 
quantum system. In Sec. IV] we study the micromotion 
effect in a trapped atom-ion system by considering the 
micromotion-induced coupling and computing the exact 
quasienergies of the system. We also suggest a simple 
scheme to bypass the micromotion effect to realize some 
interesting applications based on interaction of trapped 
atoms and ions. Finally, we offer conclusions and present 
some technical material in the appendices. 



II. ION MICROMOTION AND HARMONIC 
APPROXIMATION 

An excellent review of the classical and quantum dy- 
namics of a single trapped ion is given in Ref . [B] . Here 
we briefly summarize the key points in the motion of 
trapped ions and explain the reasoning behind the har- 
monic approximation. 



For an ion in a rf trap, the equations of motion along 
the X, y, z axes are decoupled. With uj the frequency 
of the oscillating potential and rrij the ion's mass, the 
Hamiltonian for the motion along the x axis is 
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where a and q are trap parameters which depend on rrii, 
Lu, the ion's charge Ze, and the characteristics of the trap- 
ping potential. The Hamiltonians for the motion along 
the y and z directions assume similar forms with different 
values for the parameters a and q. In a linear Paul trap, 
for instance, we have 
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which means the oscillating part of the trapping potential 
only exists in the radial directions {x, y). Typically, \q\ <C 
1 and \a\ < q^ . 

The ion's motion is a combination of a harmonic oscil- 
lation at a secular frequency wq and a micromotion that 
oscillates at the frequency w of the driving potential. The 
secular frequency is 
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thus. Wo ^ w; that is, the secular motion is much slower 
than the micromotion. These conclusions are drawn from 
the approximate solution 



x{t) « xq cos(wot) 1 + ^ cos(wt) 



(4) 



of the classical equation of motion with xq an arbitrary 
constant. Since the amplitude of the micromotion goes as 
the small parameter g, it can be seen as a jiggling motion 
around the overall path of the secular motion. When one 
averages over the short time period of the micromotion, 
the cos(a;i) term in Eq. Q vanishes and the resulting 
motion of the ion is a harmonic oscillation. This is the 
so-called harmonic approximation. 

Cook, Shankland, and Wells provide a quantum- 
mechanical derivation of the harmonic approximation in 
Ref. f7] . They write the wave function of the ion as 
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and insert this into the Schrodinger equation for the 
Hamiltonian of Eq. (II]) . It follows that the effective wave 
function ^(x,t) obeys the modified Schrodinger equation 
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where wq is the secular frequency mentioned above and 
the factor 7 is 
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Cook et al. argue that most of the fast time depen- 
dence of ^{x, t) is contained in the exponential factor 
and hence ^{x^t) may be treated as a slowly varying 
function of time. Therefore, one may take in Eq. ^ the 
time average over the short time interval 27r/a; of the 
micromotion; after that we are left with the well-known 
Schrodinger equation for a harmonic oscillator. This is 
the quantum-mechanical basis for the harmonic approx- 
imation. Equation (|6]) tells us that this approximation 
is valid only when the time-dependent terms have little 
effect on the unperturbed wave function. As we sec later, 
while this is all right for a single trapped ion, it may not 
hold when the ion is coupled to an external system. 

The Hamiltonian associated with the effective wave 
function <I>(a;,t) can be read off from Eq. (l6]), 



^eff(i) = 



2rn, 



'-TUiUjlX'^ + Ha 



.it), 



(8) 



where 



Hmmit) = —mi{'-fuJoyX'^ cos(2ujt) — jujolX, P} am{ujt) 

(9) 
accounts for the time-dependent contribution of the mi- 
cromotion. Here {X, P} is the anti commutator of the 
position and momentum operators. The form of HcS, 
where the time-independent part of the secular motion 
and the time-dependent one of the micromotion are sep- 
arated, makes it more convenient to work with when one 
wishes to study the effect of the micromotion. In the next 
section we make use of this transformation to investigate 
the trapped atom-ion system. 



III. TRAPPED ATOM-ION SYSTEMS 




FIG. 1. A trapped atom-ion system. O is the center of the 
atom trap which is cliosen as the coordinate origin; fi and 
r2 are the position vectors of the atom and ion, respectively; 
and d is the position vector of the center of the ion trap. 

The system is composed of an atom in a harmonic 
trap interacting with an ion in a rf trap as shown in 



Fig. [T] To simplify the problem we first consider the one- 
dimensional (ID) system in which the atom and ion are 
confined to moving along only one axis, say the x axis. 
The realistic three-dimensional (3D) system is discussed 
in Sec. VF below. Let us choose the origin of coordinate 
at the center of the atom trap. When the trap centers are 
separated by a distance d, the Hamiltonian of the system 
is 



H{t)=Ha+H,{t) + V, 
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The indices a and i label the atom and ion, respectively. 

The properties of the interaction potential T^nt are 
well-documented in Ref. [T]. When the distance r be- 
tween the atom and the ion is large, this potential has 
the form V-mtir) ~ — (ae^)/(2r^), where a is the polariz- 
ability of the atom. This long-range form results from the 
attractive force between the ion's charge and the atom's 
induced dipole. The characteristic length of the atom-ion 
interaction is defined by R* = \Jmae^ jK^ ^ where m is 
the reduced mass of the atom-ion system. This length 
indicates the range within which the atom-ion interac- 
tion potential is larger than the quantum kinetic energy 
f? I^mr^ . For the ^■^^Ba^ and ^^Rb system, R* is around 
5544 Bohr radii. 

The behavior of the interaction potential at short dis- 
tances is very different from its long-range form. The 
short-range potential is, in general, non-central and de- 
pends on the electronic configurations of the atom and 
ion. This spin dependence of the short-range potential is 
key to the implementation of a quantum phase gate pro- 
posed in Ref. [3j (here "spin" means the binary alterna- 
tive of two hyperfine states). Idziaszek et al. P show how 
to take into account the effect of the short-range potential 
by utilizing quantum defect theory where the short-range 
potential can be characterized by a single quantum de- 
fect parameter called the short-range phase . The basic 
idea is to replace the potential Vint with its long-range 
form while imposing a specific boundary condition on the 
wave function at a small distance r,nin <C i?*. In addi- 
tion, this distance Tmin must be sufficiently larger than 
the length scale set by the short-range potential, which 
is a few Bohr radii. 

To investigate the micromotion effect, we follow Cook 
et al. and make the transformation 



'^{xa,Xi,t)= exp 
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The resulting effective Hamiltonian for the effective wave 
function ^{xa,Xi,t) is the sum of a time- independent 
Hamiltonian and an oscillating term H^^{t) that rep- 



resents the ion's micromotion. 
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where we already replaced the interaction potential by 
its long-range form. Any calculation of the wave function 
must take into account the boundary condition dictated 
by the short-range phase. 

In the harmonic approximation, i?mm(0 i^ neglected 
and we only need to deal with a time-independent sys- 
tem. This approximate approach has been studied in 
great detail previously [H |2] . Our purpose is to investi- 
gate how the micromotion affects the unperturbed sys- 
tem; thus, we have to work with the full Hamiltonian of 



Eq. (12). We are particularly interested in how the state 



of the system evolves in an adiabatic process where the 
trap distance is changed slowly in time as such a process 
is vital for the implementation of the proposed quantum 
phase gate in Ref. |3j. 

Although the micromotion is represented by an oscil- 
lating term similar to an electromagnetic field, it is an 
intrinsic property of the system and cannot be switched 
on or off. Therefore, time-dependent perturbation theory 
which emphasizes the transitions between unperturbed 
states is not a suitable approach to the micromotion 
problem, which is the reason why we need to consider 
the Floquet formalism. 



IV. FLOQUET FORMALISM 

A. Floquet theory 

The Floquet formalism was first considered by Shirley 
[5] and has been developed in great depth for studying 
atomic and molecular multiphoton processes in intense 
laser fields [5]. Thus, we have an advanced mathematical 
framework ready in our hands to investigate the exact 
quantum dynamics of the trapped atom-ion system. In 
this section we outline the key points which are important 
for our study. 



The Hamiltonian of Eq. (12 1 satisfies the periodicity 
condition H{t + T) = H{t) with T = 2tt/lli. According to 
the Floquet theorem, the Schrodinger equation for such 
a periodic system. 
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adopts a special class of solutions called the Floquet solu- 
tions which can be expressed in terms of the quasienergy 
e and the Floquet wave function u{t) as 



*(i)=exp(--ei)u(t), 



(14) 



where u{t) is a function of both time and space coordi- 
nates and is periodic in time; that is, 



lit + T) =u{t). 



(15) 



A substitution of *(t) from Eq. (14) to Eq. (13) yields 

d 



H{t) - ih 



dt 



u{t) = eu{t). 



(16) 



Thus, the quasienergy and the Floquet wave function are 
respectively the eigenvector and eigenvalue of the opera- 
tor 



Hpit) ^ H{t)~ih 
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which is called the Floquet Hamiltonian. One observes 



from Eq. (14) that quasienergies and Floquet states are 



to a periodic system what eigenenergies and eigenstates 
are to a time-independent system. 



The solutions to the Floquet eigenvalue equation ( 16 ) 



has the following important Brillouin-zone-like struc- 
ture: If u{t) is a Floquet state with quasienergy e, then 
u{t) exp{ikujt) is also a Floquet state with quasienergy 
e -|- khuj for any integer k. These sates are physically 
equivalent because they belong to a unique wave func- 
tion, inasmuch as 



*(i) = u{t)e 
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ikujtl 
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In other words, Floquet states and quasienergies come in 
congruent classes modulo a;, each containing an infinite 
number of equivalent members. It is convenient to denote 
the Floquet states and their corresponding quasienergies 
by Un,fc and e„_fe where the index n indicates physically 
different classes and k different members in a class. Since 
for each quasienergy e there are equivalent quasiener- 
gies e -|- khiij, it is always possible to reduce an arbitrary 
quasienergy to a single zone [E — hjj/2, E + huj/2] for an 
arbitrary energy E. As there can be only one member 
from each class e„jt in a single zone, we need not worry 
about the redundancy of the physically equivalent states 
if we restrict our study of the quasienergies to one zone. 
Sambe [^ introduces the extended Hilbert space for all 
the square-integrable periodic functions with period T in 
which the scalar product is defined as 

{{u{r,t),v{f,t))) ^ (l/T) f dt{u{f,t),v{r,t)). (19) 

Jo 

This scalar product is a 4D generalization of the nor- 
mal scalar product in 3D. Time and space can be treated 
on equal footings in the extended Hilbert space. Since 
the Floquet states are constant vectors in this space, an 
advantage of working with it is that methods developed 
for time- independent quantum mechanics, such as the 
Raylcigh-Schrodinger perturbation method and the vari- 
ational principle can be readily generalized for the Flo- 
quet states. The generalization is straightforward thanks 



to the similarity between the Floquet eigenvalue equation 
( 16 1 and the time-independent Schrodinger equation. 



The final properties of the Floquet states we need to 
mention are the orthonormality 



and the completeness, 

y^ \Un,k)){{UnM = ^> 



(20) 



(21) 



which comes from the fact that any square-integrable and 
T-periodic wave function iplt) can be expanded as 



■0(i) ^^Un,k{t){{Un,k,i^)), 
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n,k 



which is just a generalized Fourier series of 'ip{t). 

Recall that without the micromotion the trapped 
atom-ion system has a time-independent Hamiltonian 
with well-documented eigenenergies and eigenstates [J 
|3]. If the micromotion is included, the system is a pe- 
riodic one that possesses well-defined quasienergies and 
Floquet states. When the strength of the micromotion 
Hmni{t) in Eq. (12) is reduced to zero, these quasiener- 
gies and Floquet states must approach the eigenenergies 
and eigenstates obtained by the harmonic approximation. 
The difference between the exact quasienergies (Floquet 
states) and the approximate eigenenergies (eigenstates) 
tells us how important the micromotion effect is. 

To be more specific, let us consider the Hamiltonian of 
the trapped atom-ion system which can be written as 



H{t)=Ho + H^^{t), 



(23) 



where Hq is the unperturbed, time- independent part with 
eigenenergies E^ and eigenstates 0„. When the micro- 
motion is taken into account these states will transform 
to some Floquet states M„(i) with quasienergies e„. An 
obvious way to quantify the micromotion effect is to cal- 
culate the energy difference 



SE = |e„ - E„ 



(24) 



and the deviation between the exact, oscillating probabil- 
ity density /9„(t) = |w„(i)p and the approximate, static 

density pn = I'/'nP- There are many ways to measure 
the difference between these probability densities. An 
example is 



5pn = 



dt I dx 



Pnix,t) - p'^Hx) 



(25) 



which is the average over time and space of the absolute 
difference between p„(x, t) and p„ (x). We see later that, 
for our particular problem, it suffices to calculate the 
energy difference; however, the oscillation of the wave 
function may have important effects in general. 



B. Floquet Hamiltonian 

In this section wc show how to obtain the Floquet 
states and quasienergies by the Floquet Hamiltonian 
method. Let us consider a system described by the 
Hamiltonian 



H{t) ^ Ho + Vcosiujt), 



(26) 



where T^ is a time-independent operator. The corre- 
sponding Floquet Hamiltonian for this system is 

d 

HF{t) = Ho-ih—+Vcos{ujt). (27) 

Suppose we know the eigenenergy En and eigenstates \n) 
of the unperturbed Hamiltonian Hq . Then the quasiener- 
gies and Floquet states of the unperturbed Floquet Hamil- 
tonian 
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are simply 



H^-" =Ho~ ih-^^ 



^nl = En + khu:, 



(28) 
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for any integer k. We may now obtain the matrix ele- 
ments of the full Floquet Hamiltonian Hp by utilizing 
the scalar product of the extended Hilbert space, 

{{u^°l\ Hp \ui°]^)) = {En + khw) <5„,™4,, 

+ -{n\V\m) (4j+i +4j-i)- 

(30) 

Diagonalizing the above infinite matrix gives the exact 
quasienergies and Floquet states of the periodic sys- 
tem. In practice one needs to truncate the Floquet 
Hamiltonian matrix by choosing n = l,...,7Ve and 
k = —Nf, .... Nf for some sufficiently large numbers Ng 
and Nf . Then, the linear size of the Floquet Hamiltonian 
matrix is iV = N^{2Nf + 1). 

Besides the numerical diagonalization of the Floquet 
Hamiltonian, we can also use perturbation methods [9' 
to find the approximate quasienergies and Floquet states 
when the oscillating potential is weak. According to 
the Rayleigh-Schrodinger perturbation method general- 
ized for Floquet states, the lowest-order corrections for 
the quasienergies and Floquet states are 
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{"i,i}7^{n,fc} 



Upon inserting V{t) = V cos ujt and e^ ■ = Em+jhuj into 



the above equations and evaluating the scalar products, 
we arrive at 



,(1) _ ikujt 
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A£',„„cos(wt) + ihuj sin(a;i) 



AEln - (f^) 



AEr, 



A£;2^-(M^' 



(32) 



where AE^m = E^ - E^- 

The fact that the energy shift is the same for all values 
of k implies that the Rayleigh-Schrodinger perturbation 
method preserves the w-modulo structure of each class 
of quasienergies. It is clear from the above equation that 
this perturbation approach is valid only when 



|(m|F|n)r«|(K„-S„ 



{tVjjf 



(33) 



which requires that no resonance exists between \n) and 
any other state |m), that is, \En — Em\ 7^ fujj. When 
there are two eigenstates close to resonance, one needs to 
use the almost-degenerate perturbation method which is 
discussed in the next section. 



C. Resonance 

To demonstrate clearly the important role of resonance 
in quasienergy structures, we first discuss a two-level sys- 
tem whose Hamiltonian has the form in Eq. ( 26 1 with 
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and 



V = h^ 
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(35) 



where the coupling constant r\ is small enough so that 
perturbation theory is applicable. Furthermore, we allow 
the frequency wq to be varied from zero to well beyond the 
frequency u) of the oscillating field. This system is exactly 
solvable but we choose perturbation methods as they are 
still useful when one goes to multilevel systems. The 
eigenenergy diagram and the corresponding quasienergy 
diagram of the unperturbed system is shown in Fig. p[a) , 
where each quasienergy level is indicated by its double 
index {n, fc}. Note that resonances appear as apparent 
crossings in the unperturbed quasienergy diagram, here 
drawn in four zones from —2a; to 2cj. 

Now we turn to the calculation of the quasienergies of 
the full system. To gain insight we first use the perturba- 
tion method to find the approximate quasienergies and 
Floquet states. Far from resonance, when 
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1/2, 



one may use the Rayleigh-Schrodinger perturbation the- 
ory to obtain 






and 
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with n = 1, 2 and n = [3 — (—1)"] /2. Thus, the energy 



shift is of the order 



Se = e„ t 
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which is very small if 77 <^ w. The small periodic correc- 
tion to the Floquet state means that this state is repre- 
sented by a wave function that undergoes small oscilla- 
tions around the mean value which is the unperturbed 
wave function. Therefore, the unperturbed wave func- 
tion is a good approximation to the Floquet state in this 
far off-resonance region. 

In the opposite extreme, where the system is close to 
resonance when ujq/uj « 1, one needs to use the almost- 
degenerate perturbation method [9j . We start with com- 



puting the Floquet Hamiltonian using Eq. (30). If one 



arranges the double index {n, k} such that n goes through 
1, 2 before each change in k, the Floquet Hamiltonian has 
this infinite block-tridiagonal structure: 
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which contains the 2x2 identity matrix / and the 2x2 zero matrix 0. 

I 



Close to resonance we have e\ ^-^ 



efl. The sub 



,(0) 



,(0) 



matrix corresponding to u\ j, 

from the full Floquet Hamiltonian and is found to be 



^ and U2 l can be read off 



H^ = 



-i,fe-i 
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(40) 



When the coupling between these two Floquet states is 
much stronger than the coupling to other states, which is 
usually true for two states in resonance, the quasienergies 
and Floquet states are approximately given by the eigen- 
value and eigenvector of the sub matrix. Therefore, at 
exact resonance where £1 {._^ — ^2 k ^ ^i ^^^ quasiener- 



gies and Floquet states are 



hr] 
ei.fc-i « eT -TT, 
2,k ^ 

'Ikujt 
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Because of the coupling, the two levels are pushed 
apart by an amount be sa hr\. Comparing this result with 



that for the far off-resonance region of Eqs. (36) and (37), 



we observe two important differences: First, the energy 
shift at resonance is much larger than its far off-resonance 
value; and second, while for the off-resonance regime the 
time-dependent part of the wave function is only a small 
oscillation around the dominant time-independent part, 
these parts possess comparable amplitudes at resonance. 
Thus, the wave function oscillates much more strongly 
when resonance happens and there is little resemblance 
between the unperturbed wave function and the real one. 
We conclude that the oscillating field has an important 
effect on the quasienergy structure and the characteris- 
tics of the wave function if there is resonance between 
any two levels of the unperturbed system. 

For the intermediate values of wq/oj, numerical diago- 
nalization of the Floquet Hamiltonian is needed to obtain 
the quasienergies. The result rapidly converges with the 
number f^ f of the Floquet modes. In our computation 
we choose N f = 20, which is verified to be sufhciently 
large. The result is shown in Fig. [2jb). The energy shift 
is indeed maximum at resonance as predicted by the per- 
turbation methods. 
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FIG. 2. (a) Eigenenergies and unperturbed quasienergies. 
The solid lines indicate the two eigenenergy levels, the dash 
ones show the zone structure of the unperturbed quasiener- 
gies. (b) A comparison of the exact quasienergies, represented 
by the dotted curves, with the unperturbed quasienergies. 
A typical avoided crossing caused by the oscillating field is 
marked by AC. 



The oscillating potential has another important effect 
when the parameter wq is varied slowly in an adiabatic 
process. As can be seen in Fig. [2Fb), there is an ap- 
parent crossing of the unperturbed quasienergy lines at 
resonance and hence the system started from the branch 
(1) will continue to move to the same state (1) in the op- 
posite end of the crossing. However, in the real situation 
when the oscillating potential is accounted for, the two 
curves are pushed apart by the coupling and we have an 
avoided crossing. Thus, the system will move from state 
(1) to state (2) in an adiabatic process. Therefore, the 
dynamics of an adiabatic process passing through a reso- 
nance changes completely when the oscillating potential 
is taken into account. 

We must stress that the adiabatic approximation for 
Floquet states is slightly different from that for the time- 
independent eigenstates. The mathematical framework 



of the Floquet adiabatic process for short laser pulses is 
described in Ref. ^H,. In Appendix [A| we modify these 
formulations a little so that it suits a more general prob- 
lem. A particularly useful result is the Landau-Zener 
formula for a Floquet adiabatic process which gives the 
approximate transition probability of a system when it 
passes through an avoided crossing. For an avoided cross- 
ing with asymptotic slope difference a and energy spac- 
ing 8t as shown in Fig.^b), the Landau-Zener transition 
probability is 



Pt 



LZ 



exp 



{belhf 



4 wq tan(a/2) 



(42) 



with wq denoting the time derivative of Wg- When a and 
(5e are known from the energy diagram, the above formula 
tells us how slow (fast) we need to change the parameter 
Wo to obtain an adiabatic (diabatic) process. 

The conclusion that the effect of a weak oscillating 
field is most important at resonances is also true for a 
multilevel system. One may use in a similar way the per- 
turbation methods to prove that the energy shift and the 
oscillation amplitude of a Floquet state u^^k are largest 
when there is a resonance between the unperturbed level 
M^^ and another level u^ fc„i- The spUt of these levels 
due to their coupling is 



(5e« |(?i|F|m)|. 



(43) 



Thus, a simple way to know whether the exact quasiener- 
gies structure differs by a great amount from the approxi- 
mate unperturbed eigenenergies is to look for resonances 
in the unperturbed energy diagram and then calculate 
the coupling strength |(n| V \m)\ between the resonating 
levels. 



V. MICROMOTION EFFECT 

We now show how we use the Floquet formalism to 
evaluate the exact quasienergies of the trapped atom-ion 
system and from there deduce the effect of the micromo- 
tion. As a starting point we assume the atom is fixed at 
the center of its trap. This corresponds to the situation 
when the atom is tightly trapped; that is, the atom trap- 
ping frequency is very large. We explain later why most 
of the conclusions about the micromotion effect in this 
simplified model also hold for a more realistic system. 

In view of the above assumption, the presence of the 
atom only results in an additional potential in the Hamil- 
tonian of the trapped ion. Therefore, the Hamiltonian of 
the system given in Eq. ( 12 ) is simplified to 



i7(t)=F0 + ffmm(t), 

- -fUJo{X, - d, PJ sin(a;t). (44) 



For the purpose of numerical calculation it is convenient 
to write the above Hainiltonian in the following dimen- 
sionless form: 
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IT 
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2 {x,/kY 
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X,-d 
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Xi — a liPi 



cos I 2 — T 

Wo 



sm 



LO 



Wo 



(45) 



where Ri = ymjmPjh? is the interaction length rede- 
fined for this problem, U — ^JhJrn^Q the harmonic oscil- 
lator length, and r = wo^ a new time variable. In short, 
all lengths are scaled by li and energies by hiOQ. 



A. Micromotion-induced coupling 

Before we compute the exact quasienergies of the sys- 
tem, it is worth examining whether resonances, as de- 
scribed in Sec 



IV C 



exist in our system. First we cal- 
culate the eigenenergies of the unperturbed Hamiltonian 
i^o- Following Idziaszek et al. [T] we compute, at zero 
trap distance d — 0, the eigenenergies i?„(0) and eigen- 
states |n(0)) of Ho{d = 0) using the renormalized Nu- 
merov method [H]. The states |?t-(0)) are then used to 
diagonalize Hq for other values of the trap distance. To 
account for the short-range potential, we need to fit the 
solution to the short-range form [1] 



^(xi) oc Xi sin 



R, 



V>6 



(46) 



for Xi -^ Ri. The short-range phase 1^93 can be related to 
the s-wave scattering length h by 
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FIG. 3. The unperturbed eigenergies of (a) the vibrational 
levels and (b) the low-lying molecular levels. The energy gaps 
between adjacent molecular levels are very large compared 
with those between the vibrational levels which appear as a 
continuum in the bottom plot. 



COt((ps) = - 



R,. 



(47) 



Once we know the scattering length b and hence the 
short-range phase, we may find a distance r-a^in satisfying 
ao *C Tniin ^ Ri such that the wave function given 
in Eq. (461 vanishes at rmm- The boundary condition 
^(^"111111) = is then used to carry out the Numerov com- 
putation. In our calculation we assume a single value for 
the odd and even short-range phases 1 . 

For our numerical calculation we consider the ^-^^Ba^ 
ion and ®^Rb atom system with Ri = 8927ao. The chosen 
trapping parameters are 7 = l/-\/2 (a = 0), wq = 27rx 100 
kHz, oj = 27r X 1.27 MHz; hence, U = 516 oo and Ri/li ~ 
17.2. The scattering length b is assumed to be 0.9i?i, from 
which we obtain r-^^^^ = 0.135 Z^. Although the secular 
frequency and the micromotion frequency stated above 
are unrealistically small, we choose these values to obtain 
plots with better clarity. We repeated our calculation 
with more realistic values of the frequencies and it did 
not result in any qualitative change in what we are going 



to describe about the micromotion effect. More of this is 
discussed in Sec. IV El 

The eigenenergy diagram of the unperturbed Hamil- 
tonian is shown in Fig. [3j At a large trap distance the 
spectrum is almost identical to that of a harmonic oscil- 
lator as it should be, [see Fig. IsFa)]. At small distances 
we have a spectrum of vibrational states [E > 0) where 
the ion is localized in the ion trap and molecular states 
{E < 0) where it is bounded to the atom. Note that the 
energy gaps between adjacent molecular levels shown in 
Fig. IsFb) are much larger than the ones between the vi- 
brational levels. The molecular states at large trap dis- 
tances correspond to the situations when we have the 
traps far apart but the atom and ion are close to each 
other. These states are, in general, not the initial states 
in which we prepare our system. The properties of the 
vibrational as well as molecular states are discussed thor- 
oughly in Refs. [3 [3]. 

In the proposed experiment for implementing the 



quantum phase gate [3], one starts with the atom and 
ion cooled to their ground states at a large trap distance 
and then adiabatically moves the ion trap closer to the 
atom trap. Thus, we are mainly interested in how the 
asymptotic ground level changes as the trap distance de- 
creases. From now on when we mention ground level we 
mean the energy curve that asymptotically coincides with 
the ground harmonic-oscillator level at a very large trap 
distance, that is, the curve marked by |0) in Fig. IsFa). 
The excited energy curves above the ground level are in- 
dicated by \n) with positive integers n, while the lower 
molecular energy curves are marked by negative integers. 
One notices that, when the trap distance is around some 
characteristic distance dc — 5li, the energy of the asymp- 
totic ground state begins to decrease rapidly. This is 
an indication that the atom-ion interaction is dominant 
within the range [0, dj- Furthermore, the energy dia- 
gram exhibits an avoided crossing where the ion changes 
from its vibrational state to the molecular state around 
dc- Since the interaction potential must be comparable 
to the trapping potential at dc, we may roughly estimate 
this characteristic distance as follows: Suppose the ion 
is half way between the atom trap and the ion trap; by 
equating the two potentials we have 



dr ~ 2R]'^iV^ 



5.2 /, 



(48) 



which is close to the value obtained numerically. While 
this good agreement may not hold for different values of 
Ri and U, the above estimation should yield a correct 
order of magnitude of the ratio ddh- 

In order to understand how the oscillating micromo- 
tion affects the energy of the ground level we look for 
resonances between this level and the excited states of 
the unperturbed system. In Eq. (44), the micromotion 



is represented by two terms with the first oscillating at 
twice the frequency of the second; hence, there are two 
types of resonance which we call the 2a; resonance and lo 
resonance. 

Since in our system uj = 12.7a;o, which is not an inte- 
gral multiple of wq, there is no resonance at a large trap 
distance (we see shortly that resonances at large trap 
distances are not important anyway). However, as the 
trap distance is reduced and the energy levels begin to 
deviate from their asymptotic values, the ground level in- 
evitably comes into resonance with a number of excited 
levels at different values of the trap distance. To find 
where the resonances appear, we plot the eigenenergies 
En together with En — oj in Fig. |4]ja) and En — 2aj in 
Fig. Bb) , which show that the ground level passes a to- 
tal number of ten w-resonances and fifteen 2a;-resonances 
with various excited levels as d goes from to 7li\ these 
resonances appear as crossings and are marked by the 
letters A„ and Bq. When the micromotion is included 
we expect the energy levels to be pushed apart by the 
micromotion-induced coupling and these apparent cross- 
ings to become avoided crossings. 

Next we need to consider the coupling strengths at the 
resonances. As mentioned in the previous section, the 




FIG. 4. A plot of (a) the eigenenergies En together with 
En — UJ shows the u) resonances and (b) the eigenenergies 
En together with En — 2a; shows the 2a; resonances. These 
resonances appear as crossings. 



shift in the quasienergy of the ground level can be in- 
ferred from the coupling strength |(0|y|n)|, where \n) 
is the excited state in resonance with the ground level. 
In Eq. (44 1, the micromotion-induced coupling is repre- 



sented by the operators 

Vi = -m,{-iujQf{X,~df, 

V2^-lUo{X,-d,Pi}. 



(49) 



Therefore, to study its effect, the quantities |(0| Vi_2 |?i)| 
need to be calculated at various trap distances. Since it 
can be shown that 



V2 \n) 



En — Eq 

ijhujQ 



Vi\n) 



(50) 



the second coupling strength is much larger than the first 
for sufficiently large n and thus has more important ef- 
fects. 

Figures [sja) and [5jb) show the collective behavior 
of many coupling strengths, obtained by numerical in- 
tegration, at four different values of the trap distance 
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d = 9li,5.5li,5.3li and d = 5.1 ^i. As can be seen in 
Fig. |3][a), the energy gap between the ground level |0) 
and the third excited level |3) is at most IIloq at rf = 0. 
Because this gap is still smaller than the micromotion 
frequency, we know that there is no resonance between 
the ground level and the first three excited levels. On 
the other hand, the energy gap between the ground level 
and the 26th excited level |26) is too large for resonance. 
Since the ground level can only come into resonance with 
those excited levels \n) with 4 < n < 25, only the cou- 
pling strengths in that range are important. We see in 
Figs. [5la) and l5lb) that when d — 9 It, the elements 
|(0| V^i,2 |"-)| are extremely small, which is due to the 
fact that the atom-ion interaction is very weak at this 
distance, and hence the unperturbed states are almost 
identical to the eigenstates of a harmonic oscillator. For 
these states we see immediately |(0ho| ^1,2 Ij^ho)! = for 
n > 2. As a result, there is no possible resonance with 
nonzero coupling strength and we observe in Sec. |IVC| 
that the energy shift caused by the micromotion in this 
case must be very small. In fact, the exact quasienergy 
spectrum of a single trapped ion was found by Glauber 
[HI to be 



taken into account. For our particular system we may 
predict that the micromotion effect begins to kick in at 
around d^^^ ~ 5.5 li. 

In the limit of a very large ratio lu/ljq, for instance 
lo/loq = 100, resonance is possible only with a large 
excitation number n. As shown in Fig. p[d), the cou- 
pling strength of Vi for large n is essentially zero at all 
trap distances and therefore can be neglected. On the 
other hand. Fig. l5Fb) shows that the coupling strength 
of V2 within the range d G [0, 5.5 Z^] only decreases slowly 
with n because it is proportional to the energy differ- 
ence [see Eq. (50l]. Even for n as large as 100 we have 
|(0| V2 1 100) I ~ O.lhujo at d = 5.3k which is not negligi- 
ble. Moreover, the coupling strength of V2 for the highly 
excited state n — 100 exhibits a sudden rise at the same 
distance dmm — 5.5 li. We verify that the sharp increase 
at dnim — 5.5 li also appears for all n G [3,100], which 
implies that the characteristic distance where the micro- 
motion effect becomes important is quite insensitive of 
the frequency of the trapping potential when the ratio 
lu/luq is in the interval [10,100], which contains most of 
the realistic values for current rf traps. 



e„ = (n+ -)hfi, 
where fi differs from the secular frequency wq by 



(51) 



(52) 



when |a| ^ q^. This number is very small, so the har- 
monic approximation works relatively well for trapped 
ions as long as the external interaction is absent or rela- 
tively weak compared with the trapping potential. 

The situation changes completely when the ion is sub- 
jected to a significant external interaction. At d = 5.5 /i, 
where the ion trap just crosses to the strong interac- 
tion region, the micromotion-induced coupling begins to 
increase rapidly and reach a remarkably large value at 
d — 5.3 li. This is due to the now strong enough inter- 
action potential which changes considerably the shape of 
the unperturbed eigenstates. To see how the coupling 
strengths change with trap distance we plot a represen- 
tative one between the ground level and the fifth excited 
level, that is, |(0| Vi_2 |5)|, in Fig.lslc). In addition to the 
sharp rise at d ~ 5.5 li, we see that the coupling due to 
Vi is indeed much weaker than the one due to V2, which 



is in agreement with Eq. ( 50 1 



Since at a small distance the coupling strength of the 
second term is comparable to the secular frequency ujq 
while the spacing of the unperturbed energy spectrum is 
of the order wq, we expect the energy diagram to change 
completely when the oscillating terms of the micromo- 
tion are included. More importantly, the fact that the 
strength of the micromotion-induced coupling increases 
sharply at a certain value of trap distance means that it 
can be used to forecast at what characteristic range dmm 
the micromotion effect becomes significant and must be 



B. Numerical calculation of exact quasienergies 

In this section we obtain, by numerical calculation, the 
quasienergy diagram, from which we then infer the im- 
portance of the micromotion effect. There are two widely 
used methods for computing the quasienergies of a peri- 
odic system: Floquet Hamiltonian and time-propagator 
methods [5]. The Floquet Hamiltonian method is simply 
the numerical diagonalization of the Floquet Hamilto- 



nian as described in Sec. IV B For our system, it can be 
expressed in terms of Hq and i/mm(i) in Eq. (44) as 



d 

Hf = Ha + H,^,^{t) - ih—. 



(53) 



We use the Floquet states of the unperturbed Floquet 
Hamiltonian at zero trap distance Hp (d = 0) , which 



are 



,(0) 



<l(0) = l^(0))e^'^"*. 



(54) 



as the basis to evaluate the matrix elements of the Flo- 
quet Hamiltonian Hp at an arbitrary trap distance 



,(0) 



(0) 



(«uo)ii?Fi<:,(o))) 



^ [ dte-''''^' {n{0)\ Hp \m{0)) e'^'^K 
T Jo 



(55) 



If one arranges the double indices {n, k} such that n runs 
through all the eigenstates before each change in the Flo- 
quet mode k, the Floquet Hamiltonian matrix has an 
infinite block-pentadiagonal structure. With / denoting 
the identity matrix and the zero matrix, the Floquet 
Hamiltonian matrix is 
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FIG. 5. Coupling strengths for (a) Vi and (b) V2 sd, d — Qh (diamonds), d — 5.5 h (circles), d = 5.3 Zi (squares), and d = 5.1 Zi 
(crosses). The coupling strengths between the ground state and (c) an intermediate excited state n = 5 and (d) a very high 
excited state n = 100 of Vi (diamonds) and V2 (circles) are also plotted as functions of the trap distance. 
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where Hq, Vi, and V2 are matrices whose elements are 

(^o)„,™ = [En{0) + Tn,ujlSl2\ 5nm " m.wg^ (n(O)l X, |m(0)) 
(l^i)n,™ = - "^.(7^0)' {nm {X, - df |m(0)) , 

(*^2)„,„ = {ih^oV [Em{Q) - E„m (yi)„^„ . 



(56) 



(57) 



Therefore, the Floquet Hamiltonian matrix can be ob- 
tained from the matrix elements of Xi and Xf at d = 0. 



Since we already obtained £'„(0) and |?i(0)) by the Nu- 
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merov method in Sec. VA the Floquet Hamiltonian ma- 



trix elements can be evaluated by numerical integration. 

In practice there are two types of truncation one needs 
to make. The first comes from the number Nf, of the 
eigenstates |n(0)) and the second comes from the number 
Nf of the Floquet mode k. To generate the basis w„ j.(0) 
we choose all the eigenstates |n(0)) in the energy range 
[—5000 Wo, 300 Wo] and k = —10, . . . , 10. In total we have 
Ne = 150 eigenstates and Nf — 21 Floquet modes; hence, 
the linear dimension of our Floquet Hamiltonian matrix 
is NeNf = 3150. To obtain a quasienergy diagram, we 
diagonalize this matrix for trap distances from to 10 li 
instep of Ad = 0.002 ^j. 

We also use the time-propagator method to calculate 
the quasienergies of the system. This method is based on 
Floquet's theorem which states that the time evolution 
operator satisfying the Schrodinger equation 



th^U{t,0)=Hit)Uit,0) 



has the periodicity property [10] 



U{t + T,0) = U{t,0)U{T,0), 



(58) 



(59) 



where U{T, 0) is a unitary operator whose eigenvalues A„ 
are related to the quasienergies e„ by 



An — 



(60) 



So, we need to obtain U{T, 0) and diagonalize it to get 
the quasienergies. First we set C/(0, 0) = / and compute 
U{T, 0) by using the propagating scheme 



U{t + At,0) =cxp 



t+At 



dt'H{t') 



U{t,0) + O{At^), 



(61) 

which is a second-order method. This simple propaga- 
tor is the lowest order of the Magnus propagator whose 
error bound is discussed in Refs. ^4^ J5,. For our cal- 
culation we use a total number of Ne = 150 eigen- 
states in the energy range [— 5000 wq, 300 wq] as a basis 
to compute the elements of the matrix H{t) and choose 
At = 10^'^r as the time step. The computation is car- 
ried out for trap distances within the interval [0, 10 h] in 
step of Ad = 0.002 Z^. In both the Floquet Hamiltonian 
and the time-propagator methods, we increased N^. up 
to 253 and repeated all the calculations to verify that we 
obtained sufficient accuracies. Here we present the re- 
sults for TVg = 150 as the quasienergy plot is clearer with 
smaller N^. 

We use both the Floquet Hamiltonian and the time- 
propagator methods to find the quasienergies and observe 
very good agreement in the results. However, while the 
Floquet Hamiltonian method is more insightful, the time- 
propagator method is found to give more accurate results 
for the quasienergies of the trapped atom-ion system for 
equal CPU time. In general, it is not convenient to use 
the Floquet Hamiltonian method when the range of the 



energy spectrum involved is considerably larger than the 
frequency of the oscillating field. The quasienergy dia- 
grams shown in Figs. [6] and [7] are obtained by the time- 
propagator method. 

As can be seen in Fig. jS) for large trap distances, the 
quasienergy of the ground Floquet state, marked by uq, 
agrees quite well with the eigenenergy obtained by the 
harmonic approximation shown in Fig. [31 The densely 
distributed curves with steep slopes are quasienergies of 
the Floquet states that correspond to highly excited un- 
perturbed eigenstates. Since all the quasienergies are 
returned in the zone [— fia;/2, fia;/2], which is an intrin- 
sic property of the Floquet Hamiltonian as well as time- 
propagator method, these highly excited quasienergy lev- 
els are projected down and they appear to cross the 
ground level. However, at large trap distances, say 
d}^ 5.7 k, the interactions between the ground level and 
those states are extremely weak and hence the avoided 
crossings, if they really exist, have very small gaps. 
Therefore, it is not hard to move the trap distance fast 
enough to diabatically cross these possible weak avoided 
crossings so that the system remains in a single Floquet 
state, which is the ground level in this case. 

As the trap distance gets smaller than a character- 
istic distance dmm — 5.7 Z^, the quasienergy begins to 
differ greatly from the eigenenergy obtained by the har- 
monic approximation. This is due to the sharp increase in 
the micromotion-induced coupling discussed in Sec. |V A[ 
Within this range, numerous avoided crossings with ob- 
servable gaps occur in the quasienergy diagram, many 
of which do not appear in the diagram obtained by the 
harmonic approximation shown in Fig. [3] The reason is 
that, at a small trap distance, the coupling of the sec- 
ond term shown in Fig. [5] is so strong that it changes the 
energy diagram completely. 

What does this mean for the implementation of the 
quantum phase gate? As we slowly change the trap 
distance passing through dnim, the large number of new 
avoided crossings of various sizes will lead to branching 
from the initial ground level to different levels [unless 
near-perfect experimental control of d{t) is carried out]. 
The branching of the wave function would result in 
a poor fidelity F of the phase gate, as we show in 
Appendix A that 1 — F oc Pe, where Pe is the small prob- 
ability of the ion being in an excited Floquet state at the 
end of the adiabatic process. Moreover, even if one can 
control an adiabatic passage of a micromotion-induced 
avoided crossing, for example the one marked by AC in 
Fig. pfb) , the duration it takes must be quite long since 
the energy gap of this avoided crossing is quite small. 
This long duration would lead to a slow gate speed 
and hence traversing the micromotion-induced avoided 
crossings is not desirable. 

The best strategy may be moving the ion trap toward 
the atom trap until it reaches the characteristics distance 
dmm and stop for a while with the hope that the ion picks 
up enough phase difference for the phase gate. However, 
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FIG. 6. (a) Quasienergies of the trapped atom-ion system plotted in the first zone [—huj/2, hij/2]. While the vibrational energy 
levels are affected by the micromotion at small trap distances, the energy curves for the low-lying molecular states stay in 
good agreement with the results obtained by the harmonic approximation, (b) An enlarged plot of the ground level Floquet 
state. The micromotion-induced avoided crossings appear and the energy curve begins to disintegrate around the characteristics 
distance dmm — 5.7 li. 



it is likely that the energy difference between the triplet 
and singlet states [5] at dnim is very small and a full phase 
of TT is not obtainable in a reasonably short time. This 
energy difference must be determined by experiments be- 
cause it depends on the short-range phase of the atom- 



ion interaction potential which is not known. In short, 
whether stopping at the minimum distance dnim is good 
enough cannot be answered by a purely theoretical cal- 
culation. 
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FIG. 7. Quasienergy levels obtained when only the first term Vi of the micromotion is included, the energy curve of the ground 
level is similar to the result obtained by the harmonic approximation except that the new avoided crossings appear at what 
was the 2ui resonances in the approximate energy diagram shown in Fig. [4] 



It is not difficult to think of an experiment to detect 
the characteristic distance dmm of the micromotion effect. 
One prepares the ion in the ground Floquet state in the 
ion trap initially placed very far from the atom trap. 
The ion trap is then slowly brought closer to the atom 
trap up to some distance dmin and back to the initial 
position. If dniin > dtnrm ouc cxpccts the ion to come 
back to its initial ground state; on the other hand, when 
rfmin < c?mm, branching would occur and only a fraction 
of the ion's wave function is in the ground state. 

As indicated in Fig. l6[a), while the vibrational en- 
ergy levels are affected by the micromotion at small trap 
distances, the energy curves for the low-lying molecular 
states, for instance the u_3 level, stay in good agree- 
ment with the result obtained by the harmonic approx- 
imation shown in Fig. IsFb) for all trap distances. Note 
that the quasienergy curve for the w_3 level goes out at 
the top and comes back at the bottom of the picture. 
The fact that the molecular Floquet states have smooth 
energy curves without any avoided crossing of observable 
size means that the micromotion does not greatly affect 
these states. The reason is that the energy gap between 
a low unperturbed molecular state and other states, as 
shown in Fig. p[b), are much larger than the micromo- 
tion frequency uj and hence the oscillating terms of the 



micromotion cannot couple the unperturbed molecular 
state to other states. In other words, the unperturbed 
molecular state does not come into resonance when the 
trap distance is changed from 10 h to 0, and as a re- 
sult the oscillating terms have little effect on the energy 
structure of these states. Consequently, the wave func- 
tion of these Floquet molecular states must consist of 
a dominant time- independent part which is just the un- 
perturbed wave function and a small time-dependent part 
that oscillates at the frequency w. As we show in the next 
section, the fact that these molecular states are almost 
unaffected by the micromotion offers a way to implement 
the quantum phase gate despite the difficulties caused by 
the micromotion that we discussed in the preceding para- 
graphs. 

To compare the effect of the two oscillating terms of 
the micromotion, we calculate the quasienergies of the 
system when only the first term Vi is included. The re- 
sult is shown in Fig. [7J the energy curve of the ground 
level is similar to the result obtained by the harmonic 
approximation except that new avoided crossings appear 
at what was the 2a; resonances in the approximate en- 
ergy diagram shown in Fig. 4rb). This structure can 
be explained by the lowest-order perturbation theory de- 
scribed in Sec. live] It confirms the observation we made 
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in Sec. |V A| that the first term only has a weak effect. 
Although they seem to be insignificant, the new avoided 
crossings lead to a completely different dynamics of the 
system in an adiabatic process. This means an oscillating 
term, even if relatively weak, must still be taken into ac- 
count if it results in resonances between the unperturbed 
energy levels. 

We see in Figs. l6Fa) and [t] that most of the dramatic 
change of the energy structure is caused by the second 
term V2 of the micromotion. This coupling is so strong 
that perturbation methods are no longer good enough 
to account for its infiuence. Because V2 contains the mo- 
mentum operator P which is related to the velocity of the 
micromotion, the intuition that the micromotion must be 
important because its classical velocity is large is correct 
after all. In the limit of very large ratio uj/ujo, we may 
neglect the first term Vi in the micromotion Hamilto- 
nian for the reasons explained in the last paragraph of 
Sec. |VA| In the same paragraph we also showed why 
the characteristic distance dnim should not change when 
the ratio uj/luq increases up to 100. Since Vi can be ig- 
nored when Lo/ujo ~ 100 , numerical computation of the 
quasienergies is considerably simpler in this limit. 

Last but not least, we must stress an important point 
regarding the numerical calculation of the quasienergies 
of a trapped atom-ion system. One may notice that the 
Floquet Hamiltonian and time-propagator method can 
be used with the original Hamiltonian given in Eq. (10 1. 



The reason we always work in the transformed picture 
is that the Cook-Shankland transformation reduces the 
magnitude of the time-dependent terms in the Hamilto- 
nian by a factor of ujo/(jJ and hence much faster numeri- 
cal convergence is achieved. In fact, given the same size 
of the Floquet Hamiltonian matrix, the accuracy of the 
quasienergies obtained when working with the original 
Hamiltonian is much poorer than that of the quasiener- 
gies obtained in the transformed picture. 



C. A scheme to bypass the micromotion effect 

We already showed that an adiabatic process using the 
ground level may not be practically possible because of 
the micromotion effect. However, the fact that the low- 
lying molecular states are almost unaffected by the micro- 
motion means that we can bypass the micromotion effect 
by making a transition from the ground level to one of 
the molecular levels at a distance larger than dnim before 
the micromotion becomes important. More specifically, 
we may prepare the ion in the ground Floquet state at a 
large trap distance and then move it adiabatically toward 
the atom trap. When the trap distance reaches some 
value di which is slightly larger than dmm, we make the 
transition from the ground level uq to a molecular level 
lying below whose wave function has a sufficient overlap 
with Uq. This transition can be realized by applying an 
electromagnetic field E{t). When the system is in this 
Floquet molecular state it will be "protected" from the 



micromotion and we may continue the adiabatic process 
up until d = 0. This modified adiabatic process can be 
used to implement the quantum controlled phase gate 
provided that the transition strengths for the triplet and 
singlet states are equal at di. 

In principle, we need to investigate the transition be- 
tween two Floquet states: Uq and a Floquet molecular 
state, say the one lying two levels below the ground level, 
U-3. Nevertheless, when d > di we know that the micro- 
motion has little effect and the dominant parts of these 
Floquet states are given by the unperturbed eigenstates 
|0) and |— 3) shown in Fig. [3J hence, it is possible to use 
the states |0) and |— 3) to study the transition between 
the exact Floquet states without encountering too much 
error. 

Suppose the one-dimensional system described in 



Eq. ( 12 ) is subjected to an electromagnetic field. 



E{t) = Eocos{ujft)ex- 



(62) 



The coupling energy resulting from the interaction of this 
field with the atom's induced dipole Pa = Aneoae/x'^ and 
the ion's charge e is 



l^rn = 



Aire^ae 






(63) 



If we choose the distance at transition to be di = 6/^, 
the frequency w/ of the electromagnetic field must satisfy 
the resonance condition 

nwf^Eo{Qli)-E^^{&k), (64) 

and a quick glance at Fig. [3lb) reveals that 

w/ w62wo = 27rx 6.2MHz. (65) 

We need to evaluate the matrix element |(0| Vcp |— 3)| to 
know whether the coupling is strong enough to make a 
transition in practice. At di the first term in the potential 
V^p is negligible compared with the second term; thus, 
we only need to compute the element |(0| Xi |~3)|. This 
gives 



|(0|X, 1-3)1 « 1.3 X 10-4 ?, 
The Rabi frequency of the transition is 



f7=^|(0|X, 



-3)1, 



(66) 



(67) 



and for an electric field strength of 2Vcm-^ we obtain 
a Rabi frequency of around 1 MHz which corresponds to 
a sufficiently strong transition. The numerical values re- 
quired for the frequency LOf and the field strength are 
both reasonable; thus, it is practically possible to carry 
out a transition from uq to u_3 at di. A transition from 
a vibrational state to a molecular state may also be real- 
ized by Raman transitions. 

Another interesting application of the scheme above is 
the creation of an atom-ion macro molecule. At the end 
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of the process when the atom and ion are in the state 
M_3 at d = 0, it is possible to make a transition to a 
much lower molecular state. After that we may turn off 
the atom trap and what is left is a bounded atom-ion 
molecule stored in the ion trap. The same procedure be- 
tween the resulting ion and another atom can be repeated 
to create a larger molecule as long as the stability con- 
dition for the ion trap is not violated. Such an atom-ion 
macromolecule may well be of interest in its own right. 



D. Excess micromotion 

In practice, there are two types of micromotion in an 
ion trap. The one we described so far is the intrin- 
sic micromotion which comes from the driving potential. 
The other type of micromotion is the excess micromotion 
which comes from imperfections in experimental setups. 
In Ref. |4j, Berkeland et al. give a detailed description 
of the excess micromotion and its unwanted effects in 
high-resolution spectroscopy of an ion in a linear Paul 
trap. Ideally, the excess micromotion must be eliminated. 
Here we investigate a realistic situation where a small 
amount of the excess micromotion is left in the motion 
of a trapped ion. 

Basically, the excess micromotion is caused by a uni- 
form static electric field -Edc or a phase difference (^ac 
between the ac potentials of the ion trap's electrodes. 
This phase shift results in an oscillating electric field 
i?ac sin (wt) , where w is the frequency of the driving po- 
tential. With Rt denoting the radius of the trap's cylin- 
drical electrodes, the amplitude i?ac is related to the 
phase difference (/Jac by [3] 



a single trapped ion is now 



Eac 



'-qiJ-Rtifa 



(68) 



Suppose i?dc and i^ac points along the x axis; the 
Schrodinger equation for the motion along the x axis of 



ih 



dt 



■^{x,t)=[ 



2mj dx^ 



+ -miU! X [a + 2gcos(a;i)] 



eEdcX — eEi^cXsm{Lut)>"^{x,t). (69) 



The term proportional to x'^ in the potential describes 
the driving potential in the ideal case when there is no 
excess micromotion. This potential leads to the intrinsic 
micromotion whose effect has been discussed in previous 
sections. The other terms which are linear in x describe 
the excess micromotion. 

To investigate the effect of the excess micromotion, 
we first follow Cook et al. [7 and ignore all the time- 
independent operators in the Schrodinger equation. The 
solution is then, obviously, 



*(a;,i) = *(x,0)exp| 



h 

eE, 



niiUJX sin(w<) 



+ ^^a;cos(wt)l|. (70) 



The phase factor here reduces to the Cook-Shankland 
phase in Eq. ([5]) when Eac vanishes. Again, one may 
argue that the main effect of the time-dependent terms 
in the potential is to introduce this oscillating phase to 
the wave function. Therefore, we are motivated to write 



the wave function in Eq. ( 69 ) as 



*(a;,t) = expj 



I 

h 

eE, 



q 2 ■ I \ 

-m.iUJx sin(a;t) 



+ — ^a;cos(a;i) \^{x,t) (71) 

with the expectation that the time-dependent part of 
^{x,t) is sufficiently small. A substitution of the above 
expression into Eq. ( 69 ) yields the following equation for 
$(x,t): 






ix,t) = { 



~\- -mujQX — eEdcX — 'mi{'-fujo) x cos(2a;i) -I- 2ihjujQ I x— — I — ) sm{ujt) 



2mi dx"^ 2 

7^0 



1 



eEacX siii(2ujt) + ih cos(ujt)— — h - — 

UJ uiiUJ ox znii 



(^) [cos(c^t)]'}<i>(:c,i). (72) 



I 

The last term gives rise to a time-dependent global phase where 
in the wave function with no physical meaning and hence 
can be discarded. Consequently, the Hamiltonian of a 
single trapped ion in the transformed picture is 



Hn = 



2to,- 



'-jTiiUjlX^ - eEdcX 



(74) 



H{t) ^ Hq + Hnunit) + H,^{t), 



is the unperturbed Hamiltonian, iifmm(^) is the Hamil- 
tonian of the intrinsic micromotion which has the same 
(73) form as in Eq. ^, and Hcx{t), the Hamiltonian of the 
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excess micromotion, is 

H^Jt) ^ -^—^eEi,cXsm(2ujt) - eE^^ cosfwt). (75) 

Lu rriiU) 

In Appendix B we demonstrate that the addition of 
the excess micromotion does not lead to any change in 
the quasienergies of a single trapped ion. In other words, 
the ion possesses a harmonic oscillator-like quasienergy 
spectrum as given in Eq. (51), where ^ is a function of 
uj,a,q and does not depend on E^c and E^c- However, 
the excess micromotion does lead to significant changes 
in the structure of the Floquet states. 

Now let us work out the effect of the excess micromo- 
tion for the one-dimensional trapped atom-ion system in 
which the atom is fixed at the center of the atom trap 
as described in Fig. [ij In this case, the presence of the 
atom results only in an additional interaction potential 
in the Hamiltonian given in Eq. ( 73 1 . Let us choose the 



location of the atom as the coordinate origin, then the 
Hamiltonian of the system is 



iI(i)=i/0 + ^mm(i)+i?cx(t) 



Ho =^ + lm,u'oiX, - df- eEdciX, - d)- 



ae 

2X} 



2m, 2 

i?mm(i)= - m,(7CJo)^(X, - dfcos{2ut) 
- ^ujo{Xi - d, Pi} sin{ujt), 

-ffox(i) = — -eEi,c{X,-d)sm{2ujt)-eEi,c — -cos{ujt). 
LU rriiU! 

(76) 

Thus, the contribution of the field E^c is to change the 
structure of the unperturbed cigenstates. Apart from an 
additive constant, the unperturbed Hamiltonian Hq can 
be expressed as 



Ho 



2mj 



1 



X, - d - 



eE, 



dc 



rriiul 



ae 

2X} 



(77) 



Thus the effect of the dc electric field amounts to shifting 
the trap distance by 



5d- 



eEdc 



(78) 



For a small value of £"^0: say, 0.01 Vm^ , with the val- 
ues of rrii and wo as given in the numerical calculation 
carried out in Sec. |V] the distance shift is Sd « 0.7 Z^. 
Here the distance shift is already significant and it can 
be even larger when either the electric field E^^ increases 
or the secular frequency ojq decreases. Thus, care must 
be exercised when such a uniform electric field exists in 
the trap. 

It is advantageous to make the change of variable 
d' = d + Sd to bring the form of iJg to that of the 
unperturbed Hamiltonian we considered in Eq. ( 44 1 . 



Then the eigenenergies and eigenstates of Ho at vari- 
ous values of d' can be readily taken from the calcu- 
lation done in Sec. V A[ With this transformation, the 



micromotion Hamiltonian can be written as the sum 
H'^^{t) + H^dt) + Hdcit), where the forms of H'^^{t) 
and H^c{i) are, respectively, those oi H^-^{t) and Hc^t) 
with d replaced with d', and H^dt) is 



ffdc(i) 



2TOj(7a;o)^M(X, - d') cos{2ujt) 
2^ujo{5d)Pi sin(wt). 



(79) 



So, by the change of variable we may now view the 
unperturbed system as unaffected while the dc field ef- 
fectively results in two additional oscillating terms in the 
Hamiltonian. The fact that a dc electric field leads to 
not only a displacement of the ion but also an additional 
micromotion is not a surprise because an off-centered ion 
is subjected to a stronger force of the driving rf field. 

To compare the magnitude of the excess micromotion 
caused by the ac field i?ac(^) and the intrinsic micromo- 
tion iJj',jjjj(i), we introduce the length scale ^ac associated 
with the force ei?ac by 



e^a 



miUJUJola 



(80) 



and write the Hamiltonians of the excess micromotion in 
the following dimensionless forms: 



Hac{t) la 



hlOn 



k 



X 



sin(2a;i) ^—^ cos(wt) 



(81) 



and 



HdS) 

hujo 



5d 



X, - d' 



1- 



k 



cos(2cjt) ■ 



IT 



sin(ajt) 



(82) 

By a comparision with the dimensionless form of i?^j,„ {t) 
given in Eq. (45 1 (with d replaced by d'), we see that the 



relative strength of the excess micromotion is described 
by the factors lac/h and Sd/li. From Eqs. (68 1 and (78) 
we obtain ^ac — RtVac- For a millimeter-sized trap with 
electrode radius i?t ~ 1 mm, a phase shift as small as 
ipac — 10~^ degrees yields lac/h — 6. Thus, the effect of 
the ac excess micromotion is comparable to that of the 
intrinsic micromotion even when the phase shift is very 
small. To avoid further complications due to the excess 
micromotion we have to keep this phase shift and the 
electric field E^c as small as possible. 

In circumstances when 5d/li and lac/h are not negligi- 
ble, the numerical methods described in Sec. |VB| can be 
readily applied to obtain the quasienergy diagram of the 
system. However, for a qualitative picture, it suffices to 
compute the couphng strength as described in Sec. |VA| 
When the excess micromotion is considered, we need to 
calculate six instead of two coupling strengths. Figure [8] 
shows the representative coupling strength |(0|V34|5)|, 
plotted against the original trap distance d, for the oper- 
ators 



V. 



7^0 



eEac{X,-d), 



Va = -eEa 



Pr 

: 1 

rriiUJ 



(83) 
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FIG. 8. A plot of the coupling strength |{0| V3,4 |5)| as a func- 
tion of trap distance. 



of the ac excess micromotion Hamiltonian i?ac(0! the 
two operators of iJdc(i) only differ from V34 by the factor 
2j6d/lac- There is a sharp rise at around dmm — 4.8 Z^, 
which differs from the corresponding value dmm — 5.3 li 
obtained in Sec. |V|A_ by 5d = Q.lli, as it should be. We 
observe that the coupling strengths, within the range of 
resonance, of all the operators in H[^^{t), Ha,c{t), and 
Hdc{t) also exhibit a sudden increase at dmm — 4.8 Z^. 

We may now describe the key behaviors of the 
quasienergy diagram. At large trap distances, that is, 
d > dmmi the interaction potential can be neglected and 



the quasienergy of the system is essentially the energy of 
a single trapped ion. Recall that the excess micromotion 
does not change the quasienergy of a single trapped ion, 
the asymptotic quasienergies of the system must still be 
given by Eq. (51). When the trap distance decreases we 



expect resonances to happen and micromotion-induced 
avoided crossings to occur in the quasienergy diagram. 
As a result, the qualitative behavior of the quasienergy 
curves shown in Fig. |6] still holds when excess micromo- 
tion of reasonable magnitudes exists in the system. The 
only difference expected is a shift of 5d in the micromo- 
tion characteristic distance dmm- 



E. A more realistic ID model 

The simplified model of a fixed atom interacting with 
a trapped ion considered in the previous sections helps us 
understand the quasienergy structure of a more realistic 
system. In fact, all of the conclusions we obtained from 
the simplified model are still qualitatively true when we 
allow the trapped atom to move. The Hamiltonian of 
this system is then given in Eq. (12 1. Let us consider 
the case when oja = wq. By switching the coordinates to 
those of the relative motion 



X^X,- Xa, (84) 

and the center-of-mass motion 

_maXa + mi{Xi~ d) 

^cm — \ , (,oOj 

ma + mi 
we may write the Hamiltonian of the system as 



H{t) = H,n^it) + H,,i{t) + Hi{t), 



Hcmit) — 
i/,el(t) = 



p2 

cm 

2M 
2m 



1 



M 



[M{-iUJQfXl^^ cos(2a;i) + 7a;o{^cm, i^cm} sin(wt)] , 



+ -mut (X - dY - 



2X^ 



m 
mi 



m{'yujof {X - df CQs{2ujt) + 7Wo{X - d, P} sin(cjf) 



Hi{t)^-2m{'^ujQfXc^^{X-d)cos{2ujt)-2'^ujQ X^^-^P + -- [X - d) P^^ sin(a;t), 



with the total mass M — 



M 
rrii and the reduced mass m — marrii / {rUa + mi 



(86) 



Without the micromotion, the Hamiltonians for the 
center-of-mass mode and relative mode are decoupled. 
It is the term Hi (t) caused by the micromotion that cou- 
ples these modes. It is clear that the Hamiltonian of the 
relative motion H^f,\{t) is, apart from some multiplica- 
tive constants, the Hamiltonian of the simplified model 
we considered in Sec. IV Al This was our motivation for 
the simplification in the first place. 



I 

The unperturbed eigenstates of the system are the 
product states 

Wa^^n) = l^cm) \n) , (87) 

and the unperturbed eigenenergies are the sum 

1^ 



En 



hujQ + En. 



where \n) and E^ are the unperturbed eigenstates and 
eigenenergies of the relative Hamiltonian as discussed in 
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Sec. V A The only difference is the replacement of the 



ion mass rrii by the reduced mass m and hence there is no 
qualitative change in the eigenenergy diagram as shown 
in Fig. [S] 

Now let us analyze how the quasienergy of the ground 
Floquet state looks like when the trap distance is varied 
in an adiabatic process. The asymptotic quasienergies 
of the system are the sum of the eigenenergies of the 
trapped atom, 



Ea = 



1 



hujf) 



(89) 



and the quasienergy of a single trapped ion given in 
Eq. (51). Furthermore, the ground Floquet state of the 



system is the product |0a) |'Uo(t)) of the ground harmonic 
oscillator state of the atom and the ground Floquet state 
of the ion. The asymptotic form of the ground Floquet 
state and its quasienergy are well approximated by the 
unperturbed eigenstates \ncm — 0,n — 0) and their eigen- 
values Eofi- 

When we decrease the trap distance, the micromotion 
will become important when the ground unperturbed 
state |0, 0) comes into resonance with an excited state 
\ncm,n), that is. 



E„ 



Eq,o = huj. 



(90) 



If the coupling strength |(0, 0| y |ncm,"-)| of the oscil- 
lating terms in H(t) is sufficiently large, micromotion- 
induced avoided crossings of observable size will occur 
along the quasienergy curve of the ground Floquet state. 
We need to compute the coupling strengths to under- 
stand the impact of the micromotion. For the two os- 
cillating terms in i?cm(i) they are, apart from the factor 
rui/M of order one, 

Mi-fiOof\{0,0\X^^\n,^,n)\ 

2 
= —huJo ( \/2(5„^„,2 + (^ne^.oj Snfi, (91) 



and 



7^0 1(0, 0| {Xcm,Pcm} |"-cm,«)| 



JQ (\/2(5„^„^,2 + (^n^^^oj Sn.O- 



(92) 



These couplings are nonvanishing only ii n — 0; that is, 
resonances must happen in the subspace of the center- 
of-mass unperturbed eigenstates for them to be of any 
importance. Moreover, as long as oj > 10 wq these reso- 
nances are only possible if ricm > 2, which implies that 
the coupling strengths described above always vanish at 
resonances. Hence, we conclude that the two oscillat- 
ing terms in ffcm(i) have negligible effect. This is not a 
surprise when we notice that Hc-m{t) is just the Hamilto- 
nian of a single trapped ion for which the micromotion is 
known to have little effect in the energy structure of the 
first few states. 



Similarly, for the two oscillating terms in H,-c\{t), the 
coupling strengths are nonvanishing only for resonances 
in the subspace of the relative-mode unperturbed eigen- 
states. Therefore, these couplings reduce to the ones we 
discussed in Sec. V A By repeating the numerical cal- 
culation shown in that section with rrii replaced by m, 
we find that these coupling strengths are significant only 

when d < 4.8Zrei with ?„! = J ^^ ~ 824 oq. 

Finally, the coupling strengths for the oscillating terms 
in Hi{t) are nonvanishing only if ricm = 1, which corre- 
sponds to the the following sideband resonance in the 
subspace of the unperturbed eigenstates of the relative 
mode 



En - Eq = h{0J -UJo) ■ 



(93) 



Furthermore, when Ticm — 1, the coupling strengths for 
the two terms in Hi(t) reduces, apart from multiplicative 
constants of order one, to those considered in Sec. |VD| 
and shown in Fig. [8] Again, a repeat of the numerical 
calculation with rrii replaced with m shows that these 
coupling strengths are significant only when d < 4.8 /rei- 
Since all the coupling strengths at possible resonances 
are either vanishing or only significant at trap distances 
smaller than a characteristics distance dmm — 4.8^rci, all 
of the conclusions we derived about the micromotion ef- 
fect in the previous sections are still true for the more 
realistic model where the atom is allowed to move. More- 
over, one could see the scheme mentioned in Sec. |VC 



still possible since the gap between the low lying unper- 
turbed molecular states of the relative Hamiltonian is 
much larger than the micromotion frequency and hence 
these states are protected from resonances and the re- 
sulting micromotion-induced avoided crossings. 

So the qualitative behavior of the micromotion effect 
is very much the same for the model considered here, 
which corresponds to loq = Wq, as it is for the simpli- 
fied model of Sec. |VA[ which corresponds to ujq <§:_ cja- 
Thus, we expect all the main features of the micromotion 
effect to occur for ujq ^ uja- These features include reso- 
nances, micromotion-induced avoided crossings, and the 
existence of the characteristic distance (imm- In the other 
regime where wq S> Wa, which means a very large trap- 
ping frequency for the ion trap, the ion must be tightly 
bounded to the center of its trap; so it is possible for 
the atom trap to come really close to the ion without 
sensing a very small amplitude of the micromotion. This 
fact can be seen easily by considering the extreme limit 
of Wo ~^ c», for which the trapped ion can be treated 
as a fixed ion at the center of its trap. In this case, we 
retrieve the simplified model described in Sec. |V A| with 
the atom and the ion interchanging their roles. It is obvi- 
ous that the micromotion does not enter the Hamiltonian 
for this model and the exact energy diagrams must look 
similar to the ones shown in Fig. [3] These observations 
imply that, for applications such as the quantum con- 
trolled phase gate which require the atom trap and the 
ion trap to be as close as possible, it is better to use a 
configuration with a large secular frequency wq. 
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The precise condition for the trapped ion to be consid- 
ered as a fixed ion, meaning micromotion would be of no 
importance, requires the value of the minimum trap dis- 
tance involved in each specific problem. This minimum 
trap distance is around d, 



2Ra la for the adiabatic 



quantum phase gate [see Eq. (48)]. At this trap distance, 
the shortest atom-ion distance is roughly dc — la — h- For 
the ion to be treated as a fixed ion, this distance must 
be much larger than the harmonic oscillator length li of 
the ion, which yields 



2-^ + ^ « 2 f ^ 

Ra Ra \ Ra 



2/3 



(94) 



After interchanging the subscripts i and a, one obtains 
the condition for the trapped atom to be treated as a 
fixed atom. 

For the numerical calculations in this section we use 
wq = w^ = 27r X 100 kHz and w = 27r x 1.27 MHz, which 
are unrealistically small for an ion trap. We repeated our 
calculation of the coupling strengths for the more reason- 
able values Wo = Wa = 27r X 1 MHz and uj — 27r x 12.7 MHz 
and observed the familiar sharp rising at djnm ~ 6.1Zrci- 
Hence, the micromotion effect, as we describe it, is still 
valid for these values of the secular frequency and the 
micromotion frequency. For typical trapping frequencies 
encountered in experiments, say oja = 27r x 100 kHz and 
Wo = 27r X 1 MHz, we have U k, 0.5 Zo- Since U does not 
differ greatly from la, the geometry of this realistic sys- 
tem is closer to that of a symmetric model with loq = uja, 
which is discussed in this section, than it is to a model 
with wq «; oja or wq » w^. 



F. Trapped atom-ion systems in 3D 

Again we consider the system illustrated in Fig. [l] but 
now the atom is stored in a 3D harmonic trap and the 
ion is confined in a linear Paul trap [1] . Let us denote the 
components of the position vectors by ra^k for the atom 
and ri^k for the ion (fc — x,y,z). The driving potential 
used to confine an ion in a 3D rf trap is 



where 



k 



Wk{t) = ^ [ofc -f 2qk cos(wi)] , 



(95) 



(96) 



and the parameters Uk and qk obey Eq. (pi). The axial 
trapping frequency of the Paul trap is uji^ = Lo^fajl and 
the radial secular frequencies are 

w-i,x = Wi,y = Wo, (97) 

where wo is given in Eq. ([3]). By taking into account that 
the trapping potential for the atom is 



^'^ = X1 9"^'»'^lfc^lfc' 



(98) 



one can show that the transformed Hamiltonian for the 
trapped atom-ion system is the sum 



of the Hamiltonian of a single trapped atom. 



(99) 



^a = E 



p2 



1 



2ma 2 ' ' 



(100) 



the transformed Hamiltonian of a single trapped ion. 



i?.w=E 



k L 



p2 , 

— — -f -m,uj^j,{n^k-d) 



-^niinv^/i 



i?mm(i) = - m,(7Wo)2(Xj - df cos(2wi) 

- 7Wo{Xi - d, Pi^a;} sin(wi) 

— mi(7a;o)^(Ki — d)^ cos(2a;t) 
-f7Wo{r,-d,P.,J, (101) 

and the interaction potential. 

At large trap distances, V\^t can be ignored so that 
the motions of the atom and ion are decoupled and the 
quasienergies of the system are simply the sum of the 
atom eigenenergies 



Ea^ 



k 



na,k 



^ ) ^a,k: 



(103) 



and the ion's quasienergies, as mentioned in Eq. (51) 



k—x.y 



1 



hiik 



1 



hwj 



(104) 



This asymptotic energy spectrum can be well described in 
the harmonic approximation. However, we know that at 
small trap distances the interaction potential will distort 
the unperturbed eigenstates and eigenenergies so that 
they come into resonance with each other. This is when 
the exact quasienergies differ greatly from the approx- 
imate eigenenergies and micromotion-induced avoided 
crossings appear in the quasienergy diagram. Idziaszek 
et al. investigated the unperturbed Hamiltonian for the 
cases of quasi- ID elongated traps and spherically sym- 
metric traps in Refs. pQ [3] and the energy spectrum 
they obtained is very similar to that shown in Fig. [3] 
Therefore, we expect micromotion-induced couplings to 
exhibit the same sudden rise at some characteristic dis- 
tance dmm as in Fig. [5] The quasinergy diagram of the 
3D system must also look similar to the one obtained in 
Fig. |6J where the energy spectrum changes completely 
and micromotion-induced avoided crossings appear at 
trap distances smaller than dnim- Thus, all the quali- 
tative statements we made for the simplified system con- 
sidered in Sec.|Vj3 should carry over to a realistic 3D sys- 
tem. The values of parameters of interest like the char- 
acteristic distance dmm cannot be computed with high 
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accuracies due to the lack of knowledge about the exact 
value of the short-range phase and must be determined 
by experiments. 

Because of the asymmetrical nature of the atom and 
Paul traps, the dynamics of the system in an adiabatic 
process depends on the direction along which we move 
the traps. From the form of the driving potential Vi{t) 
and the micromotion Hamiltonian iJmm (t) we notice that 
the micromotion occurs only in the radial directions. 
Along the axial direction the ion behaves like a har- 
monic oscillator. When the atom and ion are trapped 
very tightly in the radial direction it is likely that the 
micromotion effect is weaker if we move the trap along 
the axial direction. The ideal configuration to minimize 
the micromotion effect is the one for which the micro- 
motion frequency is as large as possible compared with 
the radial trapping frequency of the Paul trap, and for 
both traps the radial frequencies should be much larger 
than the axial trapping frequencies 17j. However, even 
then the micromotion cannot be ignored at small trap 
distances because the interaction potential Vint clearly 
couples the motion along the radial directions to the mo- 
tion along the axial direction and thus it is possible for 
the micromotion to be propagated to the motion along 
the axial direction. This transfer of the micromotion is 
negligible for relatively large trap distances but it must 
be significant at small trap distances, especially when 
the system transforms to a molecular state for which the 
atom and ion are close to each other. 

Our analysis from Sec. |V|4_ to Sec. [V^ implies that 
applications based on a Floquet adiabatic process, such 
as the creation of the atom-ion macromolecule and the 
implementation of the quantum phase gate, are realiz- 
able as long as we can carry out the following steps: (i) 
preparing the trapped ion in its Floquet ground state, 
(ii) moving the ion trap adiabatically to the atom trap 
and (iii) make a transition to a low lying molecular state 
at a trap distance di > d^^ before continuing with the 
adiabatic process until d ~ 0. Requirement (i) raises 
the question about whether the micromotion has any ef- 
fect on the cooling of a trapped ion to its ground state. 
The laser cooling process of trapped ions when the in- 
trinsic micromotion is taken into account is discussed in 
Refs. [Si [16], and the influence of the excess micromo- 
tion is studied in Ref. [1]. Even though the micromotion 
may result in unwanted heating of the ion, cooling can 
still be achieved by choosing the right value for the laser 
frequency. Therefore, the intrinsic and the excess mi- 
cromotion are not fundamental obstacles for the task of 
cooling a trapped ion to its ground Floquet state. 



VI. CONCLUSIONS 

In this work we show that while the micromotion can 
be neglected for single trapped ions it must be taken into 
account when the ion strongly interacts with an external 
system. In the context of trapped atom-ion interaction. 



the effect of the micromotion is most important when 
there is resonance between two eigenstates of the ap- 
proximate Hamiltonian. This happens when the energy 
separation is close to the micromotion frequency and the 
coupling that arises from the micromotion between the 
two states is significant. The micromotion leads to not 
only a shift in energy but also numerous new avoided 
crossings in the energy diagram. The magnitude of these 
effects can be inferred by (i) a calculation of micromotion- 
induced coupling strengths and (ii) an exact computation 
of the quasienergies of the system. Although the sec- 
ond method gives much more detail, the first method is 
much faster yet sufficiently accurate to predict the main 
features of the quasienergy diagram. Our results show 
that the approximation that the ion trap is a harmonic 
trap breaks down at trap distances smaller than a char- 
acteristic distance dnim- The quasienergy diagram for 
this region is far more complicated than the approximate 
eigenenergy diagram and thus the behavior of the system 
in an adiabatic process may change completely when the 
micromotion is included. The complication due to the 
presence of the micromotion may seriously reduce the fi- 
delity of the proposed atom-ion quantum gate. However, 
we suggest a scheme to overcome this difficulty by uti- 
lizing a transition from the initial ground Floquet state 
to a low-lying molecular state before the onset of the mi- 
cromotion effect. These transitions can also be used to 
create an atom-ion macromolecule which may indeed be 
of interest in its own right. 
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Appendix A: Adiabatic quantum phase gate in the 
Floquet picture 

Before we discuss the adiabatic quantum phase gate 
in the Floquet picture, we need to understand how a 
quantum system evolves in a Floquet adiabatic process 
[llj . We start with the Schrodinger equation 






(Al) 



where the parameter A, which is the trap distance in the 
specific case of the trapped atom-ion system, is varied 
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slowly in time. At each moment we have the instanta- 
neous Floquet states which satisfy the Floquet eigenvalue 
equation 



i?F(A,i)u„,fe(A,t) = e„,fe(A)u„,fe(A,t). 



(A2) 



To study the adiabatic evolution, one needs to employ a 
trick called the {t, t') method [18] in which we change the 
implicit time parameter t in X(t) to a new time variable t' 
and introduce a higher-dimensional wave function $(t, f) 
with the requirement that this wave function satisfies the 
"extended" Schrodinger equation 



"dtf 



ih—^{t,t') ^ HF{\{t'),t)^{t,t'). 



(A3) 



When the initial condition of ^(f = 0) is transformed to 
that of ^(t,t' — 0) in an appropriate way [TT], we have 



*(i) = *(t,t' = t). 



(A4) 



An advantage of using t he ex tended wave function ^(t, t') 
is the similarity of Eq. ( A3 ) with the Schrodinger equa- 
tion encountered in the normal adiabatic process for 
which the Hamiltonian only depends on time implicitly 
through the parameter A, and thus well-known perturba- 
tion techniques can be used. In this picture the variable 
t is treated as a fourth dimension of space and hence can 
be dropped from now on. 

We expresses the wave function 'l'(t') in terms of the 
instantaneous Floquet states, 



#(i')=^c„^fc(i')"n,fe(A(i')) 

dTe„,fe(A(T)) 



X exp 



(A5) 



and insert it into Eq. ( A3 ) to obtain the differential equa- 



tions for the coefficients Cn,k{t') which in turn can be 
solved approximately by perturbation methods. If the 
system starts in a single Floquet state, say uo,0: then 
first-order perturbation theory gives us 



c„,fe(i') = 



' dr''^^^ 



dr 



X M u„,fc(A(r)) 






uo,o(A(t)) 



xexp 



dT'[6„.,(A(r'))-eo,o(A(r'))] 



(A6) 



Recall that all the Floquet states Un,k correspond to a 
unique physical state w„ , to find the transition amplitude 
to the physical state u„ at the real time t we need to sum 
over all values of k and then replace f with t 



Cnit) = ^ CnM{t' = t) 



(A7) 



For the system to almost stay in a single Floquet state 
at any moment in time, we need |cn(i)| <S 1 for all n, 
which is the general requirement for a Floquet adiabatic 
process. 

Let us now consider the Floquet adiabatic process used 
to implement the quantum phase gate as considered in 
Sec. IV Al The adiabatic theorem states that if the initial 
motional state is a Floquet state Un{d{ti),ti) of the sys- 
tem, the final motional state at the end of the process 



4>[tf) = Un{d{tf),tf)ejip 



dre„(d(r)) 



, (A8) 



where it is understood that the quasienergy of the ground 
state is set to zero at all trap distances. Since the instan- 
taneous quasienergy e„ (d) depends on the interaction po- 
tential which in turn depends on the spin of the atom and 
the ion, we know that the phase accumulated in an adi- 
abatic process must be spin dependent [3|. 

In the spin subspace, one can create, by single qubit 
gates, a separable state 



lx(t.))=^|fc)cfc. 



(A9) 



fe=i 



with k = 00,01,10,11, and C1C4 = C2C3. Again "spin" 
means the hyperfine levels of the atom and the ion which 
are chosen as the qubit states. The Floquet adiabatic 
process turns the initial state. 



\^iU))=uoid{t,),ti)Y,\k) 



Ck, 



fc=i 



to 



\^{ts))^uoid{tf),tf)Y,\k) 



CfcC 



iBk 



(AlO) 



(All) 



fe=i 



Each of the spin terms \k) accumulates a different phase 
Ok since the phase in an adiabatic process is spin depen- 
dent. The final spin state is |x(t/)) — J2k=i l^)cfee*^'', 
and the adiabatic process is equivalent to the applica- 
tion of a quantum phase gate U — diagje'^'' } in the spin 
subspace. 

In a realistic situation there must be a small fraction of 
the motional state being transferred to excited Floquet 
levels at the end of the adiabatic process. We consider 
the simplified situation when only one dominant excited 
Floquet state u,, is occupied with a small probability pg. 
If we first ignore the spin degree of freedom and write 
the amplitude of the ground state and the excited state 
in their polar forms, the final motional state, apart from 
a global phase, is 



fc — — CO 



Htf) = VPoUo{d{tf),tf)e' 
with po + pe — 1 and pe ^ 1. 



/p^Ue{d{tf),tf)e^^ , 

(A12) 
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For the initial state given in Eq. (AlO), because of the 
spin dependence of the phases and 9' , the final total 
wave function is more complicated, 



\^{tf))=Y.\^u)\k)cue^'\ 



where 
Wk) = VPoUo{d{tf),tf) 



/PeU, 



,{d{tf),t 



f) 



(A13) 



(A14) 



with ak = 9'f. — 9k- Note that in the above we have 
ignored a very small spin-dependence in the amplitude 
i/p^ and ^/p^. 

A comparison of the realistic final wave function 
\^{tf)) and the ideal wave function |^(t/)) shows that 
when branching occurs the motional state is entangled 
to the spin state, and the adiabatic process no longer re- 
sults in a pure state in the spin subspace. The density 
matrix for the mixed state of the spin subspace can be 
obtained by taking the partial trace of the total density 
matrix |^(</)) ('J'(i/)|; the result is 






S*01 



where 



Afcj = ifjlVk) = 1 -Pe 



1 



Jio'k-aj) 



The gate fidelity is defined as [H] 
F = mm ■ 



^{X{tf)\p{tf)\x{tf)) 



(A15) 



(A16) 



(A17) 



where the minimization is done over all possible values 
of Ck satisfying the constraint J2k \'^k\'^ ~ ^ ■ 

Now we introduce the 4x4 symmetric matrix M whose 
elements are 



M, 



kj 



1 - cos(aA: - aj), 



(A18) 



and the column vector V with the components Vj — \cj\ 
It follows that 



{x{tf)\p{tf)\xitf))^l-p,V^MV, 



(A19) 



and the maximum value of V'^ MV, subjected to the con- 
straint J2k |cfeP = Ij is 



max{v'^MV} = (A^M^'^A) ^ 



(A20) 



with the constant vector A'^ — (1, 1, 1, 1). Therefore, the 
gate fidelity is 



F = yi~Pe (ATM-^A)-^ ~^~2 (^^^^"'^)" ■ 

(A21) 
There is no simple analytic form for A'^M ^A; however, 
using the fact that Mkj — ioi k = j and Mkj < 2 for 



k y^ j, one can show that V'^MV < | and derive a lower 
bound for the gate fidelity 



F> 



l-^Pe 



1 - ^Pe. 



(A22) 



In short, we see that 1 — F oc Pe- Thus, the gate fidelity 
decreases linearly with the occupancy of the excited Flo- 
quet state. 



Appendix B: Floquet state of a single trapped ion 

Here we consider a single ion confined in a rf trap with 
both the intrinsic micromotion and the excess micromo- 
tion. The Floquet state of an ion with only the intrinsic 
micromotion is given by Glauber in Ref. |T2|- Here we 
apply his elegant proof to derive the Floquet state of a 
trapped ion in the presence of the excess micromotion. 
The potential for the motion along the x axis is 

Vx{t) —-rriiUj'^ [a + 2qcosujt] X^ 
8 

- e [^dc + -Bac sin(cjf)] A, (Bl) 

and the equation of motion is an inhomogeneous differ- 
ential equation: 

X{t) + [a + 2qcos{(jjt)] —X{t) = Fq + Fi sin(wt), (B2) 

where Fq — e-Edc/'^ii, and Fi = ei?ac/"^i. 

The homogeneous equation is the Mathieu equation 
which possesses a special Floquet solution f{t) satisfying 
the initial condition 

/(O) = 1, /(O) - iiy. (B3) 

This special solution has the expression 

fit) = e»^V(0, (B4) 

where /i is called the Floquet exponent which depends 
on a, q, u! and (p{t) is a periodic function whose Fourier 
series is 



^(i)=5]C„ 



(B5) 



The coefficients C„ and the Floquet critical exponent fj, 
can be found by simple numerical procedures |21j . The 
value of v can be then computed from the relation 



1^ — fj, + uj y^ nCn ■ 



(B6) 



When |a|, Igl <C 1, the lowest-order approximation in a 
and q gives /i w wq, with ljq defined in Eq. ([3]), and [6] 



ifiit) 



I + ((7/2) cos(wt) 
1 + 9/2 ■ 



(B7) 
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If /(t) is the solution, /(t)* must also be a solution, 
and a general homogeneous solution has the form 



Xh = Afit) + Bfity 



(B8) 



When I a 1, 1 9 1 ^ 1 we have ^ <s:_ uj so the Floquet exponent 
cannot be an integer multiple of w. This means that 
f{t) and f{t)* are not periodic functions and hence no 
homogeneous solution is periodic. 



We now show that the inhomogeneous equation (B2) 
has one and only one periodic solution when |a|, j^l <C 1. 
It is easy to see that there is not more than one periodic 
solution: Assume Xp{t) is a particular periodic solution; 
any other solution x' (t) is the combination of Xp (t) and 
a homogeneous solution, that is, 



x;{t)^Xp{t) + Af(t) + Bf{tr 



(B9) 



Since /(i) and f(t)* are not periodic, it is impossible for 
x'p{t) to be periodic unless A — B ^ 0, which results 
in x'p(t) = Xp{t). Thus, we conclude that the inhomoge- 
neous equation has no more than one periodic solution. 
We prove the existence of the unique periodic solu- 
tion by explicitly constructing it by the variation-of- 
parameters method [20] • We look for a periodic solution 
of the form 

xp{t) = k{t)f{t) + k{tYf{ty (BIO) 

subject to the condition 

k{t)f{t) + k{tyf{tr^o. (Bii) 

By making use of the time-independent Wronskian, 

f{t)f{t)*-f{trf{t)^2zv, (B12) 

we obtain Xp{t) — Re{$(i)(^(t)*}, where (^(t) is the func- 
tion given in Eq. ( B5 ) and 



$(t)= £ Dr^ 



n— — oo 



with 



D^ 



2CnFo+iF,iC, 



n+l 



L-n-l , 



2l/(/i 



(B13) 



(B14) 



In the lowest-order approximation we have 



Xp(t) w dd 



1 -\ COs{Ldt) 



UJQ 



lacSm{ujt), (B15) 



where Sd and ^ac are defined in Sec. |VD] and luq is defined 
in Eq. (pi. It is obvious that Xp{t + T) = Xp{t) and so we 
have found the unique periodic solution of Eq. ( B2 ) . 



Glauber's elegant approach can be used to find the 
quasienergies and Floquet states of the trapped ion. First 
we introduce a new position operator that depends ex- 
plicitly on time, 



Xi{t)=X{t)-Xp{t), 



(B16) 



so that Xi [t) satisfies the homogeneous differential equa- 
tion. The Wronskian 



f{t)X,{t)-f{t)X,{t) 
is constant in time and thus the operator 



A^i 



2hv 



f{t)X,{t)-f{i)X,{t) 



(B17) 



(B18) 



is also a constant of motion, and it possesses constant 
eigenkets. The operator A and its adjoint are the analogs 
of the ladder operators for the quantum harmonic oscil- 
lator, and one can show that [v4,A^^] = 1. Let us define 
the states 



A|0)=0, 



Vn! 



|0) 



(B19) 



and demonstrate that they are indeed the Floquet states 
of the trapped ion. 

Although the eigenket |0) is time independent, its 
representing wave function (a;,t|0) is not. We have 
{x, t\ A |0) = 0, which yields 



— /W|- - /(i) [x - xpit)] - xp{t)f{t) \ {x, m = 0, 



(B20) 



The solution in its normalized form is 



x 



,i|0) = 



imi j{t) 



X exp 



2fi f{t) 



X - Xp{t)Y 



(B21) 



with s(x,t) = miXp{t)x/h. By applying the operator A^ 
repeatedly to the state |0) we obtain 



{x,t\n) 



1 /TOiiy\i/4 e*''^"^'*' r/(i)* 

1/2 \ 

[x-Xp{t)] 



n/2 



X H„ 



rrtiV 

WW 



xexp|^j||[x-x,(t)]4, (B22) 



where i?„ is the Hermite polynomial of order n. Upon 
inserting f{t) from Eq (B4| to the above expression, we 
get 



{x,t\n) — exp 



-i I n + - ) /ii 



Un{t), 



(B23) 
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with the periodic functions 



X exp • 



-,1/2 

[a; - Xp(t)] 






miV 






2h 



1 



i(yi(t) 



/i^(i) 



.W]' 



(B24) 



The structures of the wave functions (a;, t\n) suggest 
they are the Floquet states with quasienergies 



n + - ) hfi. 



(B25) 



Since ^ is independent of E^c and Sac, the quasienergies 
of a single trapped ion are indeed not affected by the 
excess micromotion. 

In situations where the excess micromotion vanishes, 
we have Xp(t) = and the wave functions of the Flo- 
quet states \n) reduce to the forms obtained by Glauber 
in Ref. [IT as expected. A large amount of the excess 
micromotion, which results in a large amplitude of Xp{t), 
will lead to strong oscillations in the wave functions M„(i). 

Let us now work out the kinetic energy of the trapped 
ion when i t is prepared in the state \n). Using Eqs. (B16) 
and (B18), we have 



P{t) = m,X{t) = y^ (/(i)*A + fit)A^) + m^ipit), 

(B26) 



and the kinetic energy associated with the state \n) is 

mi[ip{t)]'^ 



{n\ P{tf \n) 1 



2to, 



= - \ n 

2 



2 V 



(B27) 

The mean kinetic energy is defined as the time-averaged 
value {P{t)^)/2mi taken over one period 27r/w of the mi- 
cromotion. After inserting /(i) from Eqs. (B4) and (B5) 
to the above equation and taking the time average we 
arrive at 



1\ h 



{n\P{t)^\n) J 

2mi 2\ ' 2l V 



J2 Clip^ + riujf 



miXp{tY 



and when |a| ^ g^ ^ 1 we have 



(B28) 



{n\P{tY\n) 



1 



hujf) 



1 



2TOi V 2/ " 2 

For the ground state we have 



rriiUjl 



iSdy+llJ2 



{Pity 



Hloq 



1 



rriiUjl 



iSdf+llJ2 



(B29) 



(B30) 



2m, 2 2^ 

The cooling of a trapped ion to its ground Floquet state 

and the final value of its mean kinetic energy are studied 
in Ref. [TB] (without the excess micromotion). We see 
that the excess micromotion leads to an increase in the 
mean kinetic energy of the ground state, but this does 
not mean a worse cooling limit if we are only interested 
in the population of the ground state. 
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